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ABSTRACT 

We use a combination of N-body simulations of the hierarchical clustering of dark 
matter and semi-analytic modelling of the physics of galaxy formation to probe the 
relationship between the galaxy distribution and the mass distribution in a flat, cold 
dark matter universe with mean density fio = 0.3 (ACDM). We find that the sta- 
tistical properties of the galaxy distribution in the model, as quantified by pairwise 
velocity dispersions and clustering strength, can be quite different from those dis- 
played by the dark matter. The pairwise line-of-sight velocity dispersion of galaxies is 
sensitive to the number of galaxies present in halos of different mass. In our model, 
which is consistent with the observed galaxy number distribution, the galaxy velocity 
dispersion is ~ 40% lower than that of the dark matter and is in reasonable agree- 
ment with the values measured in the Las Campanas redshift survey by Jing et al. 
over two decades in pair separation. The origin of this offset is statistical rather than 
dynamical, and depends upon the relative efficiency of galaxy formation in dark mat- 
ter halos of different mass. Although the model galaxies and the dark matter have 
markedly different correlation functions in real space, such biases conspire to cause 
the redshift space correlation functions to be remarkably similar to each other. Thus, 
although genuinely biased relative to the dark matter on small scales, the distribution 
of galaxies as seen in redshift space, appears unbiased. The predicted redshift-space 
galaxy correlation function agrees well with observations. We find no evidence in the 
model for a dependence of clustering strength on intrinsic galaxy luminosity, unless 
extremely bright galaxies, two magnitudes brighter than L*, are considered. However, 
there are significant differences when model galaxies are selected either by morphology 
or by colour. Early type or red galaxies show a much stronger clustering amplitude 
than late type or blue galaxies, particularly on small scales, again in good agreement 
with observations. 
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1 INTRODUCTION 

Measurements of the clustering and peculiar motions of 
galaxies can, in principle, reveal the way in which galax- 
ies are related to the large-scale distribution of mass in the 
universe and thus provide strong tests of models of galaxy 
formation. In this paper we investigate three aspects of the 
galaxy distribution that are particularly sensitive to the 
galaxy-mass connection: the distribution of relative veloc- 



ities of galaxy pairs, redshift-space distortions in statistical 
measures of galaxy clustering, and the dependence of clus- 
tering strength on intrinsic galaxy properties. 

Most current discussion of large-scale structure takes 
place within the context of the cold dark matter (CDM) the- 
ory. It has long been known, however, that cold dark matter 
models in which the galaxies are assumed to have the same 
statistical distribution as the dark matter do not account for 
the relatively low rms pairwise velocity dispersion measured 
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in gala xy surveys. Indeed, the first N-body simulations of arc more strongly clustered than blue galaxies ( Carlberg et 



m gala xy surveys, inaeea, tne nrst l\-Doay simulations ot 
the stai idard CDM model ([Davis et al.| [1985|) gave values 2-3 



19!)" 



times higher that those that had been measured for galax- grini L998D. This is, of course, related to the fact that ellip- 



ies in the CfA redshift survey by Davis fc Peebles| ([L983D 



Lc Fevre et al. 



1996 



Willmer, da Costa & Pelle- 



This conflict stimulated the development of the " high peak" 
model of galaxy formation (Bardeen et al. L986), in which 



galaxies are assumed to be more strongly clustered than, 
and therefore to be biased tracers of, the dark matter dis- 



tribution. In this case, mass fluctuations at the present day 



and the associated peculiar velocities are weaker than in an 
unbiased model with the same level of galaxy clustering. 
Davis et al showed that galaxies identified with high peaks 
in the standard fio = 1 CDM cosmology could approxi- 
mately match the observed amplitude of galaxy clustering 
while having a velocity dispers ion distribution o nly sl ightly 
higher than that measured by Davis & Peebles (1983) (and 



tical galaxies are generally redder than spiral galaxies, and 
so this effect is due largely to the same physical processes 
that drive the morphological dependence of clustering. This 
colour dependence is not a s clear cut for studies of the an- 
gular correlation fun ction (Infante & Pritchet 1993, Landy 
Szalay fc Koci| 1996), where the red and blue populations 



Th 



may have different redshift distributions 

In order to interpret the significance of these obser- 
vational results, it is necessary to replace heuristic bias- 
ing schemes with physically motivated models that can re- 
alistically incorporate the role of the environment in the 
formation and evolutio n of galaxies. Semi-a n alytic mod- 



similar, in fact, to that in an Qq = 0.2 CDM model). 



els of galaxy formation (White fc Frenk 



observational determination of the pairwise veloc- 



White fc Guiderdoni 1993|, Lacey et al. 



al 



ity dispersion has been problematic. It is now recognized 
that this statistic is very sensitive to the contribution from 
galaxies in rich clusters, which were not fairly represented 



1994 



1991|, [Kauffmann 
ICoie 



1993 



et 



in early r edshift surveys ( Mo, Jing fc Borncr 



1993 



Zurek et 



al. ||l99|, |Mo. Jing fc Borner| [1997|, ^omcrvillc Primack ~ 



Noltheiiut 1997 , jomcrvillo, Davis fc Primack L997D. More 
recent estimates from larger surveys, typically sampling vol 



Somerville fc Primack 1998) follow the formation 
and evolution of galaxies ab initio, within the context of cos- 
mological theories of structure formation. These models have 
met with significant successes in describing and predicting 
observed pr operties of galaxies, such as faint galaxy num- 
ber counts ( White fc Frenk [1991 , [Kauffmann, Guiderdoni 



fc White|p"994 [Baugh, Cole fc Frenk||1996a|), propc rt 



umes oil the order of lO'^h ^Mpc^ [ are more robust and 



in much better agreement with one another (Marzke et 
al.j [199^ [ling, Mo fc Borner| |l99^, [Ratcliffe et al.||1998aD. 



19961 



a function of morphological type ( Kauffmann 
mann. Chariot fc White [l996| , [Baugh, Cole fc Fren: 



Kauff- 



19961) 



and the nature of Lyman break galaxies at high redshif t 



( Baugh et al. 



1998 



Governato et al. 



The current values are around ~ 500 km s 
as lar ge as the original measurement by 
([1983D. 



almost twice 



igjl . 



Kolatt et al 



19991) 



Davis fc Peebles 



A complementary way to explore the connection be- 
tween the distributions of dark matter and galaxies is to look 
for a dependence of clustering on galaxy properties such as 
intrinsic luminosity, morphology or colour. There have been 
many attempts to measure galaxy clustering a s a function 
of luminosity, often with conflicting conclusions ([Phillipps & 



In their simplest implementation, semi-analytic models 
can give limited information about the clustering of galaxies 
on large scales (Baugh et al, 1999) by applying the analytic 
de scription of the clust ering of dark matter halos developed 
by Mo fc White (1996). A more powerful technique recently 



Shanks|T987 , jAlimi, Valls-Gabaud fc Blanchardj [l98^, pavi 







developed consists of incorporating the semi-analytic pre- 
scription into an N-body simulation, whereby the clustering 
of galaxies ca n be determined directly and reliably down to 
small scales (Kauffmann, Nusser fc Steinmetz [l997[ [Kauff-' 



imetz JT 
[1999b| ," 



mann et 



et al.[gi 881 
tiago fc "K 

A number of recent measurements of the redshift-space cor- 
relation function of galaxies selected by absolute magnitude 
suggests a dependence of clus tering strength on l uminosity 
for gal axies brighter than L, (Benoist et al. 1996 , Tadros fc 
Efstathfo 
et al. 



1999a|, [K auffm ann et 
1999]) 



Benson et 



W e app lied such a technique 
enson et aLJ^OOCJ), 



and found that the 
semi-analytic model of Cole et alJ (1999) produced a galaxy 



199e , Willmer, da Costa fc Pellegrini 



1998, Guzzc 



1! 99l but see Loveday et al. 1995 and Hoyle et al. 1999 



for counter-examples) . 

A dependence of clustering on morphological type is, 
however, well established, with elliptical galaxies found to 
be more strongly clustered tha n spiral galaxies as mea 
sured both in angular catalogues (JDavis fc Gellei |l976|, Gio- 



correlation function in a ACDM cosmology (i.e. a CDM 
model with Q,o — 0.3 and A/3Hq = 0.7) that matches re- 
markably well the real-space correlation f unction measured 
from the APM galaxy survey (Baugh 1996). On small scales, 
the model galaxies are less clustere d than the da r k mat ter. 
A similar conclusion was reached by Pearce et al. (1999) for 
galaxies formed in c qsmolo gical gas dynamics simulations 
and by Colin et al, (1999) for galactic dark matter halos 



vanelli, Haynes fc Chincarinj |198(: j_ [ovino ct al 



in three- dimensional surveys ( jMoorc ct al. 1994, Loveday et 
al. 199! i). This result can be understood as a reflection of 



993) and 



the morphology-density relation, the statement that ellip- 
tical galaxies preferentially inhabit high density regions of 
the Universe (Dressier 198C| ). The dependence of clustering 
on the colour of galaxies has also been measured in red- 
shift surveys with the general consensus that red galaxies 



We write the Hubble constant as Hq = WOh km s ^ Mpc ^ . 



formed in high-resolution N-body simulations. 

In this paper, we combine our semi-analytic model of 
galaxy formation with high-resolution N-body simulations 
to investigate the clustering properties of galaxies selected 
according to various criteria. We consider clustering in real 
and redshift-space and show that galaxy velocity dispersions 
are significantly lower than those of the dark matter, in rea- 
sonable agreement with the data. Our results for this partic- 
ula r statistic differ so n iewhat from those recently obtained 
by Kauffmann et al, ( 1999a ) and we explore the reasons 
for these differences. Kauffmann et al. also investigated the 
dependence of clustering amplitude on galaxy properties, 
finding that both early type and red galaxies were more 
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strongly clustered than the ga laxy population as a whole 



in a greemen t with the work of 
metz 



KaufTmann, Nusser & Stein- 



( '- 997 ) and with the trends seen in the observational 



data. However, they found no evidence in their models for 
the luminosity dependent clustering a mplitude measured by 



Willmer, da Costa & Pellegrini (1996). Our r esults for these 



pr operties are quite consistent with those of Kauffmann ct 
al. (1993a), but we compare our models to different datasets. 
The rest of the paper is laid out as follows. In §0, we 
briefly describe the semi-analytic model of galaxy formation 
and the N-body simulation used in this study. In ^ we de- 
scribe how we take into account the effects of redshift space 
distortions and calculate galaxy velocity dispersions. We ex- 
amine the model correlation functions of galaxies selected 
by luminosity, morphology and colour and compare these 
results with the available observational data in M. Finally, 
in jjo w e present our conclusions. 



DESCRIPTION OF THE NUMERICAL 
TECHNIQUE 



We implement the semi-analytic model of Cole et al. ( 199!: ) 
in the "GIF" N-b ody simulations , a full descri ption of which 
ma y be fo und in Jenkins et al, (1998) and ^auffmann et 
|al.| (1999a). These are high-resolution dark matter simula- 
tions that follow 17 million particles within a cosmological 
cubical volume of side 141.3/i~^ Mpc. In this paper, we fo- 
cus attention on the ACDM simulation, which has cosmo- 
logical parameters, fio ~ 0.3, A/{3Hq) — 0.7, h = 0.7, a 
power spectrum shape parameter F = 0.21, and normalisa- 
tion as = 0.9 (to match the local abundance of hot X-ray 
clusters; White, Efstathiou & Frenk 1993, Eke, Cole & Frenk 
1996.) Dark matte r halos are identifi ed using the friends-of- 
friends algorithm ( Davis et al. 1985 ) with a linking length 
of 0.2 times the mean interparticle separation. We consider 
only halos consisting of ten or more particles, which Kauff- 



mann ct al. (1999a) have shown are stable clusters in the 
GIF simulations, and correspond to a lower halo mass limit 



of 1.4 X lO"/i"^M0. Benson et al, (2000) have shown that 



galaxy catalogues produced by the semi-analytic model with 
this halo mass limit are complete for galaxies brighter than 
Mb — 51og/i « —18.3. The semi-analytic model populates 
the simulated volume with approximately 12,000 galaxies 
brighter than A/b — 5 log h = —19.5. 

In order to investigate the clustering properties of very 
bright, low abundance galaxies we have used a different sim- 
ulation of a much larger volume. This was carried out by 
the VIRGO Consortium and hereafter we will refer to it 
as the "512"^" simulation. It has identical cosmological pa- 
rameters to the GIF ACDM simulation, although the initial 
conditions were generat ed using the transfer function cal- 
culated from CMBFAST ( Seljak fc Zaldarriaga 1996|) rather 
than the fitting function of Bond & Efstahiou (1984) which 



was used for the GIF simulations. The simulated volume is 
a cube of side 479/i~^Mpc, populated with 134 million par- 
ticles, giving a particle mass of 6.8 x lO^'^h'^MQ. Resolving 
halos of 10 or more particles in this simulation results in 
galaxy catalogues which are complete for galaxies brighter 
than Mb — 5 log h « —20.3. Thus, we will only use this sim- 
ulation for studying the distribution of the very brightest 
galaxies. 



We have found that the clustering strength of the dark 
matter in the 512^ simulation is slightly larger than in 
the GIF simulation at large pair separations. The ratio 
of the real-space dark matter correlation functions in the 
two simulations, C5123/CGIF, is 1.4 at pair separations of 
10/i~^Mpc, decreasing to 1.1 at 4/i~^Mpc, and becoming 
unity at 2/i~^Mpc. This difference is due primarily to finite 
volume effects in the GIF simulation, which is missing some 
of the large-scale power included in the 512'^ simulation. 
This deficit is carried through into the galaxy correlation 
functions. It should be noted, ho wever, that thi s effec t does 



not invalidate the conclusions of Benson et al. 



(200C) and. 



if anything, it slightly improves the agreement between our 
models and the APM correlation function on scales greater 
^10/i"^Mpc. 

We use the same implementation of the se mi-analytic 
galaxy fo rmation model described in detail by Benson et 
al. ( 200C| ) to populate our simulations with galaxies. Briefly, 
for each dark matter halo identified in the simulation at the 
redshift of interest, a Monte-Carlo merger history is gener- 
ated, b ased upon th e extended Press & Schechter formal- 



ism (Press fc Schechtcr 



1974j^|Bondet al. 



1991, 



Bowoi 



1991 



Lacey & Cok 1993| , Cole et al. 1999). A set of simple, phys- 
ically motivated rules that describe the galaxy formation 
process is then applied to the merger history. The physics 
modelled include: 

(i) the radiative cooling of virialised gas in halos; 

(ii) the depletion of the resulting reservoir of cold gas by 
star formation; 

(iii) the feedback arising from stellar winds and supernovae 
that regulate the supply of cold gas; 
(iv) the chemical enrichment of the gas and stars; 

(v) the mergers of galaxies due to dynamical friction. 

The star formation history of each galaxy is then used to 
determine its magnitude in any required band. The effect of 
dust e xtinct ion is determined using the models of 



Ferrara 



ct alj ( |1999| ) as described by |Cole et alj (|999|) where full 
details of the semi-analytic modelling can be found. 

Our fiducial model is essentially the same as that of 



Cole et al. (1999), whose parameters are set by the require- 



ment that the model should reproduce, as closely as possible, 
a variety of properties of the local galaxy population, with 
particular emphasis placed on the shape and amplitude of 
the local, bj-band lumi nosity function mea sured in the ESO 
Slice Project (ESP) by |Zucca et al.| ( |l997| ). Two cosmolog- 
ical parameters, a s and F, d if fer sl ightly in this work from 
the values used by |Cole et al. (1999). This is simply because 
Cole et al. had the freedom to choose these values (they se- 
lected as ~ 0.93 and F — 0.19) whereas we are restricted 
to using the values employed in the GIF ACDM simula- 
tion. These small differences have no si gnificant ef f ect on 
any of the results presented here. Since Cole et al. (199!:) 
used the Press-Schechter mass function for dark matter ha- 
los, whereas we use the mass function produced by the GIF 
ACDM simulation, we altered slightly the star formation ef- 
ficiency, e*, in the model in order to obtain as good a match 
to the luminosity funct ion as Cole et al. did. Specifically, we 
set e* = 0.01, whereas |Cole et al.| ( |l999| ) chose e* = 0.0067 
The semi-analytic model identifies the most massive 
galaxy in each halo as the central galaxy, onto which radia- 
tively cooling gas is focussed. This galaxy is placed at the 
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centre of mass of the simulated halo, and is given the pecu- 
liar velocity of the centre of mass. Satellite galaxies in the 
halo are assigned the position and peculiar velocity of ran- 
domly chosen dark matter particles within the halo. In this 
way, satellite galaxies always trace the dark matter within 
a given halo. 

The technique that we em ploy differs in two impor - 
tant respects from that used by Kauffmann et al. (1999a). 



Firstly, Kauffmann et al. extract the full merger history of 
each dark matter halo directly from the N-body simulation. 
Although this approach has the advantage that the evolu- 
tion of any individual galaxy can be followed through the 
simulation, the properties of the merger histories generated 
by the Monte-Carlo method are statistically equivalent to 
those obtained directly from the N-body simulation (but 
have superior mass resolution). In particular, the merger 
histories in the N-body simulation are found not to corre- 
late with environment, in ag reement with the assumptions of 
the Mo nte-Carlo ap proach (Somerville et al. 1998, 



Lemson 



& Kaulfmann 1999). As a result, any statistical properties 
such as the correlation function, are not affected by the dif- 
ferences between the two approaches. Secondly, Kauffmaim 
et al. adjusted the parameters of their semi-analytic model 
to match the zero-point of the TuUy-Fisher relation between 
spiral galaxy luminosity and rotation speed rather than the 
observed luminosity function as we do. As a result, some of 
their models g ive a poor match to the luminosity function. 
(2000) have argued that a good match to the 



Benson et al 



luminosity function is required for a robust estimate of the 
correlation function which, furthermore, for ACDM provides 
a good match to the data. 

Our semi-analytic model also follows the evolution of 



galaxy morp hology, as d e scrib ed in Baugh, Cole & Frenk 
(1996t) and Cole et al. (1999), a feature that we exploit 
in j^.2 to investigate the dependence of clustering strength 
on morphology. In the model, virialised gas cools to form 
a rotationally supported disk and quiescent star formation 
turns this gas into stars. A galactic bulge can be built up 
through the accretion of relatively small satellite galaxies or 
as the result of a violent merger between galaxies of com- 
parable mass. In the latter case, any material in a disk is 
rearranged into a bulge and a burst of star formation occurs 
if cold gas is present. The morphology of a galaxy can thus 
change with time and is related to its environment, which 
determines both the rate at which gas cools and the merger 
rate between galaxies. 



3 REDSHIFT SPACE DISTORTIONS 

The clustering of galaxies in real space cannot be measured 
directly unless a distance indicator is available that does 
not rely on redshift. However, the real space clustering can 
be inferred indirectly, by deprojecting either the correlation 
function measured in terms of projected galaxy separa tion 
or the angular correlation function (e.g. [Peebles 198C). In 
the case of noisy projected measurements, both methods 
require assumptions about the form of the two-point corre- 
lation function in real space, whilst the second method also 
requires knowledge of the redshift distribution of galaxies. 

The three dimensional correlation function in redshift- 
space can be measured directly, by inferring the radial posi- 



M 1 - 




Figure 1. Real and redshift-space correlation functions in the 
ACDM model. The redshift-space correlation function of model 
galaxies brighter than Mb — 51ogh = —19.5 is shown as a solid 
line, while the redshift-space correlation function of dark matter 
is shown as a dotted line. The dashed and dot-dashed lines show 
the real space correlation functions for galaxies and dark matter 
respectively. Error bars on the galaxy correlation functions are 
derived from Poisson statistics. The symbols with error bars are 
observational determinations of the galaxy correlation function: 
triangles show the redshift-space correlation function of galaxies 
brighter than Al-g - - 51ogh = —19.5 from the ESO Slice Project 
(Guzzo et al 1999), whilst squares show t he real-spac e galaxy 
correlation function from the APM survey (Baugh L99q). 



tion of each galaxy from its redshift, assuming a uniform 
Hubble flow. However, galaxies have peculiar motions in 
addition to the Hubble flow, induced by inhomogeneities 
in the local density field. This leads to a modification of 
the apparent pattern of galaxy clustering, commonly re- 
ferred to as "redshift-space distortions". Generally, these 
distortions reduce the strength of clustering on small scales, 
where the motion of galaxies within virialised structures 
tends to smear out galaxy positions inferred from redshifts 
(creating 'Fingers of God'), and boost the strength of clus- 
tering on large scales where galaxies experience coherent 
flows into ove r dense regions an d out of underdense regions 



( Kaisei 



1987 , Hamilton 



1992). The redshift-space galaxy 



correlation function in our model therefore depends on both 
the spatial distribution of the galaxies and their distribution 
of peculiar velocities. We will comment on uncertainties that 
this introduces in 



i.2 



The redshift-space clustering of galaxies in the simu- 
lation is straightforward to compute from the position and 
peculiar velocity of each galaxy. Redshift space positions 
are calculated by adding v^/Ho to the x coordinate of each 
galaxy, where v^ is the x-component of the peculiar veloc- 
ity of that galaxy. This gives the redshift space position (at 
z = 0), as viewed by an observer located at infinity. From 
the resulting catalogue we can compute clustering statistics 
in the usual way. 

In Fig. hi we show the two-point correlation function 
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of galaxies brighter than Mb — 5\ogh = —19.5 in redshift 
space. The figure also shows the real-space correlation func- 
tion for comparison. The two effects of redshift-space dis- 
tortion described above are readily apparent in this figure. 
In real space the galaxy correlation function has an almost 
power-law form between 0.1 and 10h~^ Mpc, with a slope of 
7 ~ —1.8. In redshift-space, the galaxy correlation function 
falls below this power-law on small scales, whereas on scales 
larger than a few megaparsecs, it rises above it. 

Kaiser (1987) showed that the amplification in redshift- 



space of the two-point correlation function on large scales 
due to coherent inflow onto large-scale structures is approx- 
imately a factor of 1 -I- |/3 -I- i/?^, where (3 = 0,o'^/b and b 
is the galaxy bias, assumed to be independent of scale. For 
the model considered here, Slo ~ 0.3, and the galaxies are 
essentially unbiased on large scales. The amplification fac- 
tor is therefore ~ 1.37. Whilst Kaiser's expression assumes 
Ao = 0, in practice a non-zero cosmological constant makes 
very little di fference, merely in creasing the amplification fac- 
tor to 1.39 (Lahav et al, 1991). The value measured for our 
galaxies at separations of 10h~^ Mpc, where linear theory 
is still a fair approximation, is 1.27 ± 0.04, in reasonable 
agreement with the analytic expectation. 

It is remarkable that whilst galaxies and dark matter 
have very different two-point correlation functions in real 
space, on scales smaller than 3h~^ Mpc the redshift-space 
correlation functions of galaxies and dark matter are almost 
indistinguishable. This implies that galaxies and dark mat- 
ter have different velocity dispersions on these scales, and 
suggests that estimates of galaxy bias obtained from clus- 
tering in redshift space should be treated with caution. 



3.1 Galaxy peculiar motions 

Pairwise velocity statistics for dark matter and galaxies are 
plotted in Fig. gI The data points in the figure were obtained 
from the Las Campa nas Redshift Survey (h e reafte r LCRS, 
Shectman 6^11 [l996|) by |ling. Mo fc Bornei] (|l99^ ) . FoUow- 
ing [Jenkins et aL[ (1998), a quantity that may be compared 
to these data is the projected, line-of-sight velocity disper- 
sion given by: 



'^los('"p) 



J ar)4roi{r)dl 



(1) 



jar)dl ' 

where rp denotes projected separation, I distance along the 
line of sight, r = ^/rp + P, and the integral is taken along 
the line-of-sight between ±cxd, although in practice conver- 
gence is attained with integration limits of ±25/i~^ Mpc. 
The quantity o-^^^j is the line-of-sight pairwise dispersion, 
which is given by: 



rl{vl)/2 + imvf,)-{v,,f 



rt, + P 



(2) 



,1/2 



where {y\j ' " is the rms relative pairwise velocity perpen- 
dicular to the vector connecting each pair, ri2 , («n ) is the 
rms relative pairwise peculiar velocity along ri2 and (i'2i) is 
the mean relative pairwise peculiar velocity along ri2, all of 
which are functions of r = |ri2|. 

The line-of-sight velocity dispersion of the dark mat- 
ter, shown in the left-hand panel of Fig. bl is significantly 



larger than the value for LCRS galaxies on small scales. 
At r ~ 0.5/i~^Mpc, the observed value is 535 ± lOOkms"^, 
whereas for the dark matter (Tioa ~ 870kms^'^. For dark mat- 
ter, we find that the measured crios from the GIF simulation 



is in good agreement with that obtained by Jenkins et al 



( 1998) from a simulation of a volume approximately 5 times 



larger than the GIF simulation and from that in the 512 
simulation, differing by less than 50kms~^ on all scales con- 
sidered here. 

To compare to the LGRS data we have constructed a 
sample of galaxies brighter than Mr — 5 log h — —21.5 which 
corresponds approximately to the apparent magnitude limit 
of the LCRS at its mean depth of 300h~^ Mpc. The line-of- 
sight velocity dispersion of these galaxies, shown in the right- 
hand panel by the thick solid line, is in reasonable agreement 
with the data. A sample selected to have Mb — 5 log h < 
— 19.5 has a very similar crios, as shown by the thin solid line 
in the right-hand panel. 

Fig. g shows that model galaxies display a markedly 
different velocity dispersion from the dark matter on all the 
scales considered. This "velocity bias" is statistical in origin, 
and is unrelated to dynamical friction in clusters, because 
each galaxy in our model is assigned the peculiar velocity 
of a particle belonging to a dark matter halo in the N-body 
simulation. This effect arises because the galaxy distribution 
does not constitute a Poisson sampling of the dark mat- 
ter distribution. Firstly, L* galaxies are found only in dark 
matter halos more massive than ~ 10^^h~^ Mq, whereas a 
significant fraction of the dark matter is in smaller mass ob- 
jects. Secondly, the number of galaxies in each halo is not 
in direct proportion to the halo mass, as w ould be required 
for Poisson sampling (Benson et al. 2000). The number of 



galaxies per halo determines how well the velocity structure 
of the halo is traced by the galaxy distribution. Our results 
for galax y velocity dispersions di ffer somewhat from those 
fou nd b y Kauffmann et al. ( 1999a ) using a similar technique. 
In §3.2 we explore thoroughly the causes of these differences. 
The kind of biases seen in the pairwise velocity disper- 
sion function are also present in single galaxy statistics. Fig. 
ra shows the three-dimensional velocity dispersion of galaxies 
and dark matter within spheres of radius r, measured rela- 
tive to the mean velocity within the sphere. At large r, this 
quantity tends to a constant value, the rms peculiar veloc- 
ity of individual galaxies. As the figure shows, the dispersion 
within spheres is lower for galaxies than for dark matter and, 
furthermore, it depends on galaxy luminosity. For bright 
galaxies, the velocity dispersion is about 20% lower than 
for the dark matter. The dependence on luminosity arises 
because the fainter sample contains relatively more galaxies 
in the most massive halos which have the largest velocity dis- 
persions. The figure also shows that the sampling variance 
of the velocity dispersion function is very large, particularly 
on small scales. For example, for sphere radii of 10h~^ Mpc, 
velocities as low as ~ 100 km s~^, are within 1-a of the mean 
dispersion. This reflects the importance of long wavelength 
density fluctuations in determining the small-scale velocity 
held. 



3.2 Comparison to previous work 

The pairwise galaxy velocity dispersions that we find are 
somewhat lower than those obtained by Kauffmann et al 
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Figure 2. Pairwise velocities of dark matter particles (left-hand panel) and galaxies (right-hand panel) measured in the ACDM simulation. 
The dotted line is the mean inward radial peculiar velocity of pairs, < i;2i >; the short-dashed line is the rms pairwise radial peculiar 



velocity, < v'^. >^'^; the long-dashed line is the rms pairwise perpendicular peculiar velocity, < v'^ >^^^; the solid lines are the rms 
pairwise line-of-sight peculiar velocity, crioa; and the dot-dash line is the Hubble expansion given by iiuub = —Hof, where Hq is Hubble's 
constant and r is physical separation. In the right hand panel the thick solid line shows trios for a sample with Mr — 5 log h < —21.5. 
whilst all other lines are for a sample of galaxies with Mb — 51og/i < —19.5. The data points with error bars are taken from Jing, Mo 
& B6rn([r| ([ ' 



3) and show the pairwise velocity dispersion, (712, estimated for the Las Campanas redshift survey. The error bars are the 
square root of the sum in quadrature of the errors derived from the data and Itr uncertainties inferred from mock catalogues by Jing et 
al. These points should be compared to the line-of-sight dispersions, fTi^g, for the models (solid lines). Data points and ai^^ are plotted 
against projected separation, rp, whilst all other curves are plotted against true separation, r. 



(1999a) using exactly the same N-body simulations as us, 
but a different semi-analytic prescription. Kauffmann et al 
found very similar (to within 10%) pairwise velocities for 
dark matter particles and galaxies on a ll scales. Thus, un- 



l ike ou rs, their model disagrees with the Jing, Mo fc Bornci 
( 1998 ) data. Since the two semi-analytic treatments attempt 
to model the same physics, it is instructive to understand 
the reasons for this discrepancy. 

The first, and most obvious difference is that Kauff- 



mann ot al. (1999a) considered galaxies one magnitude 
fainter than we do and the pairwise velocities are higher 
for fainter galaxies (because fainter samples consist of pro- 
portionally more satellite galaxies in clusters than brighter 
cuts). In our model, the line-of-sight velocity dispersion, ctios, 
of galaxies brighter than Mb — 5 log ft — —18.5 is approxi- 
mately 150 km s~^ higher than that for galaxies one mag- 
nitude brighter. However, as may be seen in Fig. Ml, even for 
galaxies of the same luminosity, there is still a discrepancy 
between our results and those of Kauffmann et al. There 
are three effects that could give rise to this discrepancy: (i) 
different assignments of peculiar velocities and/or positions 
to galaxies in the same halo; (ii) differences in the frequency 
with which the two semi-analytic models populate a halo of 
given mass with a particular number of galaxies; and (iii) 
statistical variations in the formation histories of galaxies in 
halos of a given mass. We have carried out several tests of 
these effects. 

The results of the tests are shown in Fig. M. The thick 



solid line shows aios for the dark matter in the GIF sim- 
ulation, whilst the thin solid line shows aios for galaxies 
brighter than AIb — 5 log ft < —19.5 in our model. In or- 
der to estimate the variance arising from different galaxy 
formation histories, we have generated ten realisations of 
our model using different random number sequences when 
constructing dark matter merger trees and selecting random 
halo particles on which to place the galaxies. The variation 
in (Tioa between these ten models is an underestimate of the 
true statistical uncertainty because the population of dark 
matter halos is, of course, identical in all ten realizations. 
The dashed line shows crios for galaxies in our model brighter 
than Mb — 5 log ft < —18.5, with la error bars derived in this 
way. Although aios for these galaxies is larger than for our 
brighter sample, it is still biased relative to the dark matter 
and is significantly below the Kauffmann et al. result. Thus, 
statistical uncertainties and magnitude selections, although 
relatively large, do not, on their own, explain the difference 
between the two semi-analytic models. 

We now consider differences in the way in which the two 
models assign positions and velocities to galaxies within a 
given halo. In our model, each galaxy is assigned the position 
and velocity of a randomly chosen dark matter halo parti- 
cle (except the central galaxy which is assigned the position 
an d velocity of the c e ntre of mass of the halo). In the model 



of Kauffmann et al. (1999a), galaxies initially form on the 
most bound particle in their halo. If that halo later merges 
with a larger one, the galaxy becomes a satellite in the new 
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Figure 3. The three-dimensional velocity dispersion of individual 
galaxies and dark matter in spheres of radius r measured relative 
to the mean velocity within the sphere. Filled squares show the 
velocity dispersion of the dark matter, whilst open symbols show 
the velocity dispersion of galaxies brighter than Mb — 5 log h = 
— 19.5 (triangles) and A/b — 51ogh = —18.5 (circles). The symbols 
indicate the mean velocity dispersion, and the error bars show the 
rms scatter in this quantity. 



halo arid remain s attached to t he da rk matter particle on 

(1999) have shown that this 



Diaferio et al 



which it formed, 
placement scheme leads to a dynamical velocity bias, such 
that blue galaxies in clusters have a higher velocity disper- 
sion than the cluster dark matter because they preferentially 
occupy particles which are in the process of falling into the 
cluster for the first time. Such dynamical biases cannot (by 
construction) arise in our model. 

To test the effect of the different placement schemes, we 
have taken the galaxy catalogue of Kauffmann et al. (lim- 
ited at Mb — 51og^ = —18.5 and kindly provided by A. 
Diaferio) and identified the dark matter halo of the GIF 
simulation to which each of these galaxies belongs. We have 
then assigned positions and velocities to these galaxies us- 
ing our own placement scheme. The resulting ctios is shown 
in Fig. W as the dot-dashed line (with error bars calculated 
from ten realisations of the placement scheme). This may 
be compared to the crios for the original Kauffmann et al. 
catalogue shown as the dotted line. The two curves are very 
similar on small scales and the small differences that there 
are on scales larger than lh~^ Mpc are, in fact, due to our 
choice of always siting one galaxy at the centre of mass of 
each h alo. (When we place all the Kauffmann et al. galaxies 
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Figure 4. Line-of-sight velocity dispersions, ai^s for: dark mat- 
ter (thick solid line) ; galaxies in our model that are brighter than 
AfB - 5 log h = -19.5 (thin solid) a nd Mr - 51ogfo = -18.5 
(dashed line); galaxies in the model of Kauffmann ct al. (1999a) 
that are brighter than Mb — 51og?t = —18.5 (dotted line); the 
Kauffmann et al. galaxies but with positions and velocities as- 
signed using our placement scheme (dot-dashed line). The error 
bars on these lines reflect sampling uncertainties in the simulated 
volume as discussed in the text. Squares with error bars are es- 
timated from the Las Campanas Redshift Survey (Jing, Mo & 
Bornei||l99^ . 



ticular number of galaxies in the two semi-analytic models. 
This statistic is shown in Fig. g. Our model and that of 
Kauffmann et al. agree well for halos below ~ lO^^/i~^M0, 
but for higher masses Kauffmann et al. assign systemati- 
cally more galaxies to each halo (typically around 30% more) 
than we do. Thus, the catalogue of Kauffmann et al. gives 
greater weight to the most massive, highest velocity disper- 
sion halos to which pairwise statistics are quite sensitive. 
As a result, they predict significantly higher velocity disper- 
sions than we do. Since the IQ^^ h~^ Mq bin contains only 
three halos in the GIF simulation, we also show our model 
predictions (filled symbols) from an independent sample of 
twenty halos of mass IQ^^ h~^ Mq populated using exactly 
the same galaxy formation rules (indicated by filled sym- 
bols). For comparison, we have indicated in the figure the 
number of galaxies brighter than Mb — 5 log ft = —19.5 
found in the Coma cluster and the corresponding number 
in our model. The virial mass and radius of the Coma clus- 
ter are 0.8 ±0.1 x 10"/i.~^A/Q and 1.5/i.~ ^Mpc (1.2°) re- 
spectively (Seller, Diaferio & Kurtz 1999). Within a cir- 



cular aperture of this radius centred on NGC4874, God 



on ranc omly chosen halo particles the two curves agree ex- 



tremely well on all scales.) This test demonstrates that the 
differences between our results and those of Kauffmann et 
al. are not due to the different ways in which galaxies are 
identified with dark matter halo particles in the two models. 
We conclude that the differences between our results 
and those of Kauffmann et al. must be due to the frequency 
with which halos of different mass are populated with a par- 



win, Metcalfe & Peach (1983) find 100 galaxies brighter 



than A^B — 5 log ft = —18.5 and 25 galaxies brighter than 
A/b — 5 log ft = —19.5 (at the fainter magnitude cut a con- 
tamination of approximately 5.5 field galaxies is expected). 
Although this sample of one cluster cannot be used to con- 
firm our model, it is consistent with our predicted distri- 
bution of galaxy numbers but also with that predicted by 
Kauffmann et al. To differentiate between these two models 
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Figure 5. The number of galaxies per halo brighter than M-g — 
5 log ft = —1 8.5 as a function o f halo n iass in our model (squares) 
and that of Kauffmann et al. (1999a) (triangles). Also shown is 
the number of galaxies per halo brighter than Mg — 5 log h = 
— 19.5 as a function of halo mass in our model (circles). Filled 
symbols show the corresponding predictions from our model for 
an independent sample of twenty 10^^h~^MQ halos. The symbols 
show the median of the distribution whilst error bars indicate the 
10% and 90% intervals of the distribution. The number of galaxies 
brighter than Mb — 5 log ft = —18.5 a nd —19.5 within the viria' 
radius of the Coma cluster as found by [Godwin, Metcalfe fc Peach 
(1983) are indicated by stars. 



would require a larger observational sample of rich clusters. 
Our model has been shown to be in good agreement with 
the observed total luminosities of gr o ups a nd clusters de- 
termined by [Moore, Frenk fc White] ( |l993| ) (Figure 16 of 
200(]D, indicating that it is populating both 



Benson et al 



groups and clusters with galaxies in a realistic fashion. 

Although velocity statistics are sensitive to the numbers 
of galaxies placed in high mass halos, we have checked that 
galaxy correlation functions (in both real and redshift-space) 
are much less affected. For example, the redshift-space cor- 
relation function of galaxies brighter than Mb — 5 log h < 
— 18.5 differs by only 10% on scales of 10/i~^Mpc and by 
less than 20% on scales of 0.3/i~^Mpc when estimated using 
the numbe r of galaxies per halo pred icted by our model and 
by that of Kauffmann et al. (1999a). This is comparable to 



th e scatter between mo dels with different parameters found 
by [Benson et alJ (pOOCl). 



Fig. o shows the effects of varying key model parame- 
ters on the pairwise velocity dispersions. In 



4.1 



we describe 

these models further and and explain our criteria for re- 
jecting models as being non-realistic. Of the fifteen variants 
considered, five would be classed as non-realistic, and are 
shown as dashed lines. These models show a greater spread 
in velocity dispersions than the viable models (solid lines) 
which have only a small scatter around the fiducial model. 
This demonstrates that robust predictions for velocity dis- 
persions can be made once non-realistic models are excluded. 
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Figure 6. The line-of-sight pairwise velocity dispersion, ai^s, for 
galaxies with Mg — 5 log ft < —19.5. The lines show fifteen vari- 
ants of our fiducial model in which some of the key parameters 
have been varied. The five models shown as dash ed li nes would 

The data 



be rejected as being non-realistic, as discussed in S4.1. ^ 

points with error bars are taken from [ling. Mo fc Borner (1998) 



and show the pairwise velocity dispersion, 0-12, estimated for the 
Las Campanas redshift survey. 



4 THE DEPENDENCE OF CLUSTERING ON 
GALAXY PROPERTIES 

We now examine the clustering properties of various sub- 
samples extracted from our original galaxy catalogue, se- 
lecting by luminosity, morphology or colour. 



4.1 Dependence of clustering on luminosity 

In Fig. M we show our predicted redshift-space correlation 
functions for four cuts in absolute magnitude, and compare 
them with the corresponding correlation functio ns measured 
in the ESP redshift survey by Guzzo et al. (1999). (The space 
density of galaxies in our model at each of these magnitude 
cuts matches closely that found in the ESP survey, since 
our bj-band luminosity function is tuned to agree with the 
ESP survey — c.f. Fig[ll[). The error bars on the model are 
estimates of the sample variance, computed in the following 
manner. For each magnitude cut, we extract fifty randomly 
placed cubes from our simulation which have a volume equal 
to that of the corresponding volume-limited samples drawn 
from the ESP survey. For the ESP survey volume-limited 
at Afbj — 5 log ft, < —18.5, —19.5 and —20.5, our simulation 
volume contains 46, 17 and 6 independent survey volumes 
respectively. Thus, for the brighter samples our estimates 
of the sample variance is likely to be an underestimate of 
the true value. In the figure, we plot the median correlation 
function of these fifty realisations, and error bars showing 
the 10% and 90% intervals of the distribution. As can be 
seen, the model results are in excellent agreement with the 
data over the whole range of scales and magnitudes shown. 
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Figure 7. Redshift-space correlation functions of galaxies se- 
lected by bj-band absolute magnitude (see panels for magnitude 
cuts). The lines show the median model correlation function ob- 
tained from fifty randomly chosen volumes equal in size to the 
corresponding ESP sample. The error bars on these lines indi- 
cate the 10% and 90% intervals of the distribution of correlation 
functions. The model result for Mb, — 51og/i < —19.5 is repro- 
duced, for reference, in each panel as the dot-dashed line. The 
filled squares with error bars show the correlation fun ctions cs- 
ti mated from the volume-limited samples of the ESP ( puzzo et| 
al. ||l999l . 



Guzzo et al. (1999) argue that for the brightest mag- 
nitude cut they considered, Mbj — 51og/i < —20.5, there 
is evidence for stronger clustering than for galaxies with 
Mbj — 5 log h < —19.5. Unfortunately, this effect is most no- 
ticeable on scales somewhat larger than those we can probe 
reliably with our simulation. At r ~ 10/i~^Mpc, the largest 
scale on which we can reliably measure the correlation func- 
tion, we find that galaxies with Mb, — 51og/i < —20.5 have 
a clustering amplitude that is only 1.07 ±0.10 times greater 
than that of galaxies with Mbj — 51og/i < —19.5 (when 
measured from the full volume of our simulation). Thus, the 
results from the GIF simulation are inconclusive. We have 
therefore made use of the 512"^ simulation, described in §S, to 
examine the clustering of very bright galaxies. The correla- 
tion function of galaxies brighter than Mbj —5 log /i = —20.5 
in this simulation agrees with that of the same galaxies in 
the GIF simulation, within the rather large errorbars (due 
to the limited number of such bright galaxies in the GIF 




Figure 8. Redshift-space correlation functions of galaxies 
brighter than Afb — 5 log h = —20.5 (solid line) and -21.5 (dashed 
line) measured in the 512^ simulation. The redshift-space corre- 
lation function of dark matter in this simulation is shown by the 
dotted line. 



volume), except on the largest scales where the different 
dark matter correlation functions in the simulations (as dis- 
cussed in ^2) introduce a similar difference in the galaxy 
correlation functions. As shown in Fig. pi we find that, in 
the 512"^ simulation, galaxies with Mbj — 5 log ft < —21.5 
have a clustering amplitude on scales above 10h~^ Mpc 
which is approximately 4 times higher than that of galax- 
ies with A/b, — 5 log ft < —20.5. The space densities of 
Mbj — 5 log ft < —20.5 and -21.5 galaxies are approximately 
6 and 200 times lower than that of L» galaxies respectively. 
The large increase in clustering amplitude is due to the fact 
that these very bright galaxies are found mostly in halos of 
mass ^Af», for which the halo bias is a rapidly changing 
function of halo mass (Mo & White 199(:). 

In Fig. y we show the real-space correlation function 
for o ur models app l ying t he same four absolute magnitude 



cuts. Guzzo et al. (199£) estimated the real-space corre- 
lation function by fitting a power-law to the correlation 
function measured in terms of projected separation. We 
plot their results over the range of scales used in their fit, 
0.5 < r < 9.5ft~^Mpc. Within the estimated errors from 
sample variance, our models are in reasonable agreement 
with the measurements. Guzzo et al. claim to find evidence 
for a weak luminosity dependence in the real-space corre- 
lation function, with galaxies in the brightest magnitude 
cut being the most strongly clustered. This effect is man- 
ifest as an increase in the correlation length with intrinsic 
luminosity obtained from their power-law fits to the correla- 
tion function. The slope of the ESP correlation function also 
increases for intrinsically brighter galaxies. In contrast, we 
find no significant evidence for a luminosity dependence of 
either clustering length or the slope of the correlation func- 
tion in our models. The pair-counting errors that Guzzo et 
al. estimated for the fitted parameters are largest for the 
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Figure 9. Real-space correlation functions of galaxies selected 
by bj-band absolute magnitude (see panels for magnitude cuts). 
The dashed lines show the power-law fits t o the correspondin g 
ESP survey correlation functions obtained by puzzo et al. ( 199!^ ). 
Dotted lines track the quoted errors on these fits. The solid line in 
each panel shows the median model correlation function obtained 
by splitting our sample into fifty randomly positioned cubes with 
volume equal to that of the ESP survey at the same magnitude 
cut. Error bars on this line indicate the 10% and 90% intervals of 
the distribution from these same samples. The model result for 
Mbj — 5 log ft < —19.5 is reproduced, for reference, in each panel 
as the dot-dashed line. 



brighte ;t sample, which contains the fewest galaxies. Our 



model correlation functions, however, indicate that sample 
variance is particularly significant for the fainter catalogues, 
which sample smaller volumes. Taking this sample variance 
into account weakens the significance of the effect claimed 
by Guzzo et al. 

To confirm that our mode l pred ictions for ^(r) are ro- 
bust, we follow Benson et al. (2000) and consider the ef- 
fect on the correlation function of galaxies of varying key 
parameters in our semi-analytic model, for example, the 
galaxy merger rate or the baryon fraction. Each model is 
constrained to match the ESP bj-band luminosity function 
at L* by adjusting the parameter T, which gives the ratio of 
the total mass in sta rs (including brown dwarfs) to that in 



luminous stars. (See Benson et al, (EOOOl) for a description 



of these variant models. 

In Fig. hd we show the real-space correlation functions 



Figure 10. Real-space correlation functions of galaxies selected 
by bj-band absolute magnitude (see panels for magnitude cuts). 
The fifteen variant models are plotted either as dotted lines if they 
produce a correlation function similar to the standard model or 
as dashed lines if they do not. 



obtained from the fifteen models with altered parameter val- 
ues. Of the fifteen variants considered, five exhibit signifi- 
cantly different clustering properties to the fiducial model. 
As can be seen in Fig. hll, these models are amongst the 
ones that disagree the most with the observed luminosity 
function, where these models are shown as da shed lines. 
This hig hlights again the point emphasized by (Benson et 



al 



2000) that a good fit to the galaxy luminosity function 



is a pre-requisite for a robust determination of clustering 
statistics in semi-analytic models. The same five models were 
shown as dashed lines in Fig. H, were they show a greater 
variation in velocity dispersions than models which are a 
good fit to the luminosity function. This shows that repro- 
ducing the bright end of the luminosity function is also im- 
portant in order to make robust predictions for galaxy ve- 
locity statistics. 

4.2 Dependence of clustering on morphology 

In Fig. [13, we compare the real-space correlation functions 
for galaxies selected by morp hological type in our m odel 



with observational data from Loveday et al, (1995). Al- 
though we restrict attention to galaxies in the same mag- 
nitude range used by Loveday et al., our criteria for mor- 
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Figure 11. The local bj-band luminosity function. Symbols with 
error bars show various observational estimates of this quantity, 
as indicated in the legend. The lines show fifteen variant models 
along with our standard model (solid line). The five models which 
show significantly different correlation functions are identified by 
dashed lines. 
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Figure 13. The mass function of host dark matter halos, 
weighted by the number of galaxies with —20.0 < AI^, — 51og/i < 
— 19.0 per halo, for late type galaxies (solid line) and early type 
galaxies (dotted line) in our model. The mass functions are nor- 
malised by the total number of galaxies in each sample. The ar- 
rows indicate the mean host halo mass for each population of 
galaxies. 
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Figure 12. Real-space two-point correlation function for galaxy 
samples selected by morphological type. Squares and triangles 
with error bars show the observed correlation functions of early 
and late typ e galaxies respe c tively , with —20.0 < Mb, — 5 log h < 



19.0, from Loveday et al. (1995). Solid and dashed lines show 



correlat ion functions for model galaxies in the same magnitude 



range. "! 'his sample is subdivided by morphological type as de- 
scribed m the text, lor the model, the lines indicate the corre- 
lation function from the whole simulation volume, which is ap- 
proximately 1.3 times larger than the M^^ — blogh < —19.0 
volume-limited Stromlo-APM survey. The dotted line shows the 
correlation function of the dark matter. 



phological classification are inevitably somewhat different. 
Loveday et al. classified galaxies by eye as either early types 
(E/SO) or late types (S/frr). Our classification is based on 
the bulge-to-total luminosity ratio of each galaxy as mea- 
sured in dust-extincted B-band light (B/T). We classify 
galaxies with B/T< 0.4 as late type galaxies and those with 
B/T> 0.4 as early type galaxies. The observed correlation 
between bulge-to-total luminosity ratio and morphological 
type display s a large scatter (see Fig. 1 of Baugh, Cole & 



Frenk 



1996b) 



Bearing these caveats in mind, the level of agreement 
between the observed correlation functions and the model 
predictions shown in Fig. [13 is encouraging. Early type 
galaxies are the most stron gly clustered on a ll the scales 
we consider, as found also by 



Kauffmann et al. 



( [I999a| ). The 



greatest difference occurs on scales below Ih^ Mpc, where 
early type galaxies have a clustering amplitude that is more 
than 4 times higher than that of late type galaxies. On larger 
scales, the difference in clustering amplitudes persist, but is 
less pronounced. 

In Fig. nsl we plot the host halo mass functions of late 
and early type galaxies in our model (weighted by the num- 
ber of galaxies per halo) . As may be seen, early type galaxies 
do, on average, reside preferentially in higher m ass environ- 
ments th an late type galaxies, as speculated by Loveday et 
al. (1995). This explains the larger correlation length of the 



early types sinc e higher mass ha l os are intrinsically more 
strongly biased ( |Frenk et al] [1984 |Mo & Whit^|l99e| ) rela- 



tive to the dark matter distribution as a whole. It also ex- 
plains the greater correlation amplitude of early type galax- 
ies on scales below lh~^ Mpc, where the correlation function 
is sensitive to the number of galaxy pairs inside cluster sized 
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Figure 14. The real-space correlation function of galaxies 
brighter than Mb — 51og/i = —19.5 in our model selected by 
B-R colour. The solid line is the correlation function of galaxies 
with B-R> 1.3, whilst the dashed line is that of galaxies with B- 
R< 1.3. Error bars on these lines show the Poisson pair counting 
errors. The dotted line shows the correlation function of the dark 
matter, and the points with error bars are the real-space corre- 
l ation functio n measured from the APM survey for all galaxies 
( Baugh||l996| ). 



halos (Benson et al. 200C|) . Late type galaxies are rare in 



these halos and so their pair count on these scales is small, 
resulting in a lower correlation amplitude. 

4.3 Dependence of clustering on colour 

At present, the only available observational constraints on 
the colour dependence of galaxy clustering come from deep, 
magnitude-limited samples that cover a significant baseline 
in redshift. The technique used in this paper needs to be 
extended to include evolution both in the galaxy population 
and in the clustering of dark matter to produce realistic 
mock catalogues that can be compared with these observa- 
tions. We defer a detailed discussion of these techniques and 



a thorough comparison to the data to a later paper (Benson 



et al 



2000 ). However J forthcoming redshift surveys, such a s 



the 2dF (CoUesE 1996) and SDSS (Gunn fc Weinberg 1995) 



will contain large enough numbers of galaxies with photom- 
etry in different bands to allow the construction of local 
volume-limited samples subdivided by galaxy colour. 

We present, in Fig. fl4l, our predictions for colour de- 
pendent clustering of galaxies in real space, measured in the 
full volume of our simulation. Red galaxies are defined to 
have B-R> 1.3 and blue galaxies B-R< 1.3, wh i ch m atches 
the cut used by [Willmer, da Costa fc Pellegrini| ( |l998|) . The 



error bars show the Poisson pair counting errors in each bin. 
The clustering of red galaxies is significantly stronger than 
that of blue galaxies, which follows from the fact that the 



red galaxies are predominantly early types, and so are pre- 
dicted to reside preferentially within higher mass halos than 
the blue galaxies. 



5 DISCUSSION 

In this paper, we have combined N-body simulations of the 
hierarchical clustering of dark matter with semi-analytic 
modelling of the physics of galaxy formation, to probe the 
relationship between the distribution of galaxies and the dis- 
tribution of mass in a A-dominated cold dark matter uni- 



In an earlier paper (Benson et al. EOOOl), we studied the 
correlation function of galaxies brighter than L* in real-space 
(i.e. as a function of true spatial separation; see also Kauff- 
mann et al 1999a). We found remarkably good agreement 
between the predictions of the ACDM model and measure- 
ments of the correlation function of the APM galaxy survey 
(Baugh 199C), over four orders of magnitude in correlation 
amplitude. On small scales, we found the model galaxies to 
be less strongly clustered than (or biased low relative to) the 
dark matter. We have now seen that the net effect of small- 
scale peculiar velocities is to cancel out this bias: in redshift 
space the galaxy and dark matter correlation functions are 
very similar to one another. This cancellation arises because 
the pairwise velocity dispersions of galaxies and dark matter 
are different. Fortuitously, the difference is just sufficient to 
compensate for the differences in real-space clustering. Thus, 
although genuinely biased, the distribution of galaxies mea- 
sured in redshift space appears, to a good approximation, 
unbiased on small scales. The pairwise velocity dispersion 
of the model galaxies is ~ (200 — 300) km/s lower than 
that of the dark matter over almost two decades in spatial 
separation and is in g ood agreement with rece nt observa- 
tional determinations ( Jing, Mo fc Bornei 1998). However, 
compar ison of our results with those of |Kauffmann et al 



( 1999a ) demonstrates the sensitivity of this statistic to the 
number of galaxies that populate rich clusters. Both models 
are consistent with the available data for the Coma cluster, 
but produce line-of-sight velocity dispersions which differ by 
around 100 km/s. 

The physical origin of the offset between the galaxy and 
dark matter velocity distributions lies in the way in which 
galaxies sample the velocity field of the dark matter. The 
mass-to-light ratio of halos in our model is a strong, non- 
monotonic function of halo mass. Galaxy formation is most 
efficient in halos of mass ~ 10 h~ Mq, and the mass-to- 
light ratio increases at higher and lower masses due to long 
cooling times f or the gas and stron ger feedback respectively 
(see Fig. 8 of Benson et al. MOO). As a result, the num- 
ber of galaxies per halo does not increase as rapidly as the 
halo mass. Thus, when computing the pairwise velocity dis- 
persion, high velocity dispersion halos are undersampled by 
galaxies relative to the contribution of these halos to the 
velocity dispersion of the dark matter itself. 

We have explored the sensitivity of our theoretical pre- 
dictions for peculiar vel ocities to variations in the galaxy 

( ^000| ) we found that 



Benson et al. 



formation parameters. In 
predictions for the galaxy two-point correlation function in 
a particular cosmological model are robust to changes in 
model parameters, provided that the galaxy luminosity func- 
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tion re mains approximately the same. Here, we find a similar per and also in KaufFmann et al.| ([1999aj) and [Benson et 



result far the peculiar velocities; models with similar lumi- 



nosity functions produce similar results. The amplitude of 
the velocity bias depends not only on the luminosity func- 
tion, but also on the shape of the power spectrum of density 
fluctuations and on the statistics of the occupation o f ha- 

(pOOOl), our 



Benson et al. 



los by galaxies. As we showed in 
(cluster-normalized) ACDM model gives a good match to 
the observed luminosity function of galactic systems. To a 
large extent, this is the reason why the model also gives 
a good match to the observed pairwise velocity dispersion 
function (although our neglect of dynamical biases in the 
galaxy distribution may affect our results slightly) . The dif- 
ferences bet ween our predicted pecul iar velocities and those 
obtained by Kauffmann et al. (1999a) from the same N-body 
simulations, but using a different semi-analytic model, are 
simply due to differences in the way in which the two models 
populate high mass halos. These, in turn, reflect differences 
in the luminosity functions of galaxies and galactic systems 
in the two models. 

The dependence of clustering on intrinsic galaxy prop- 
erties provides, in principle, an interesting test of models 
of galaxy formation. The most obvious property to consider 
is galaxy luminosity. Unfortunately, the dependence of clus- 
tering on luminosity is difficult to measure from magnitude 
limited redshift surveys and so the observational situation is 
inc onclusive. A weak effect has been claimed, for example, 
by Guzzo et al. ( 1999 ) for galaxies approximately one magni- 
tude brighter than L» . For the bulk of the galaxy population, 
we find no evidence for a dependence of clustering on the 
intrinsic luminosity of our model galaxies. This is not sur- 
prising in view of the fact that galaxies of a given luminosity 
reside in dark matter halos spanning an appreciable range of 
masses. We do, however, predict a strong effect for galaxies 
that are significantly brighter than L,. Unfortunately, the 
space density of these galaxies is too low for their clustering 
to be measured in existing redshift surveys. On the other 
hand, both models and observations agree that early type 
or red galaxies are more strongly clustered than late type 
or blue galaxies, particularly on small scales. This is a re- 
flection of the observed morphology-density relation which 



ar ises naturally in hierarchical clustering models ( Frenk et 



al. 



198^ , [Kauffmann||l996| , [Baugh, Cole fc Frenk| p^gebU 
ft is worth emphasizing that the galaxy formation 
model that we have used in this p aper is ess e ntiall y the 



same as that discussed at length in Cole et_al 



ce pt for the smaU differ ences mentioned in 
by [Benson et al| ([200(]|). 



( [1999| ) (ex- 
and adopted 



Cole et al. flxed parameter values 
by requiring their model to reproduce a variety of proper- 
ties of the local galaxy population, with emphasis placed on 
achieving a good match to the local bj-band galaxy lumi- 
nosity function. This same ACDM model also reproduces 
the present day clustering of galaxies in real and redshift 
space, is in reasonable agreement with the observed line-of- 
sight pairwise velocity dispersion of galaxies, and matches 
t he clustering meas ured at intermediate and high redshifts 
( [Baugh et aij[l999| ). 

An important conclusion of this work is that the sta- 
tistical properties of the galaxy distribution can be quite 
different from those of the underlying dark matter. A phys- 
ical approach to modelling the formation and evolution of 
galaxies, using the kind of techniques discussed in this pa- 



al. (2000), is therefore imperative if we are to understand 
the biases in the way that different galaxies trace the dark 
matter. This, in turn, is a pre-requisite for making sense of 
the unprecedented amount of information about the galaxy 
distribution in the low and high redshift Universe that will 
shortly become available from forthcoming redshift surveys. 
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